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We consider the condensation energy, Ec, in superconductors where the pairing is electronic in ori- 
gin and is mediated by a collective bosonic mode. We use magnetically-mediated superconductivity 
as an example, and show that for spin-fermion couplings A ~ 1, the physics is qualitatively different 
from the BCS theory as the condensation energy results from the feedback on spin excitations, while 
the electronic contribution to Ec is positive due to an "undressing" feedback on the fermions. The 
same feedback effect accounts for the gain of the the kinetic energy at strong couplings. 
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Undertanding the origin of the condensation energy is 
an important step towards identifying the mechanism of 
high temperature superconductivity in the cuprates. In 
a BCS superconductor, the condensation energy - the 
energy gained upon entering the superconducting state, 
is E^'-'^ — ~NfA'^/2, where A is the superconducting 
gap and Nf is the fermionic density of states. This 
decrease in the total energy upon pairing results from 
a fine competition between an increased kinetic energy 
and a decreased potential energy, both of which are much 
larger than Ec- For optimally doped YBCO with Nf ~ 
10~'^TOeT^~"'^, and A ~ ZQmeV , the BCS formula yields 
Ec ^ 5K in reasonable agreement with Ec ~ 0.12fcs7c ^ 
lOK that Loram et al [gj extracted from specific heat 
measurements. However, at smaller dopings, A increases 
while the measured Ec decreases, in clear disagreement 
with the BCS theory. 

In this paper, we present a computation of the con- 
densation energy for a strongly coupled superconductor. 
We argue that at strong coupling, the relation between 
Ec and A is different from BCS theory, and is consistent 
with experimental trends in the underdoped cuprates. 
We show that the agreement with the BSC theory at 
optimal doping is largely coincidental as at this doping 
range, strong coupling effects are already dominant. 

Some earlier works have already suggested that the 
gain in the condensation energy in the cuprates may be 
due to a non BCS physics. Scalapino and White Q con- 
jectured that at strong coupling, the dominant contribu- 
tion to the condensation energy comes from a feedback 
effect on the magnetic excitations of the system. On 
the other hand, Norman et al argued Q that the con- 
densation energy likely has an electronic origin, and is 
driven by a gain in the kinetic energy which at strong 
coupling is negative (in contrast to BCS theory) because 
of a strong "undressing" of fermions which bear a greater 
resemblance to free particles in the superconducting state 
than they do in the normal state. Similar ideas were 
expressed by Hirsh and Marsiglio ||^. A different idea, 
related to the lowering of the Coulomb energy in the su- 
perconducting state has been proposed by Leggett Q . 

In this communication we argue that these apparently 
disparate viewpoints are in fact consistent with each 



other, and describe the same strong coupling physics. 
In the superconducting state, spin decay into fermions is 
forbidden at energies smaller than 2A. This simultane- 
ously gives rise to two effects. First, the spin propagator 
develops an excitonic (resonance) peak at LUres < 2 A. 
The energy released by the creation of an exciton re- 
sults in a gain in the magnetic part of the condensation 
energy. Secondly, the fermions cannot decay until their 
frequency exceeds A -|- ujres (this is the magnetic ana- 
log of the Holstein effect). The elimination of fermionic 
scattering at low frequencies implies that the fermionic 
self-energy T,(lu) in the superconducting state is reduced 
compared to that in the normal state. This effect low- 
ers the kinetic energy (see below) and at strong coupling 
overcomes the effect of particle-hole mixing that increases 
Ekin in a BCS superconductor. 

The above reasoning is quite general and can be ap- 
plied to various superconducting mechanisms. A more 
subtle issue is to explain the reduction of Ec in the un- 
derdoped regime, and why Ec lOK at optimal dop- 
ing despite the fact that by all accounts, strong coupling 
effects are already large. Below we present an explicit 
computation of Ec assuming that the pairing is due to 
spin-fluctuation exchange We show that our theoretical 
Ec agrees with the data both in magnitude and doping 
dependence. We view this agreement as support for the 
spin-fluctuation scenario for the cuprates. 

The condensation energy Ec is the difference between 
the grand free energies in the superconducting and nor- 
mal states ||] i?c = f^s — ^n- We obtained Ec by mod- 
ifying the Eliashberg formula for phonon superconduc- 
tors 1^, ^ to the spin case. This approach implies that 
Ec is dominated by momenta near the hot spots where 
the momentum dependence of the d— wave gap can be ne- 
glected except for the sign change between "hot" regions. 
We will discuss the validity of the Eliashberg theory for 
the spin case after we present the results. In the Eliash- 
berg theory, Ec can be represented as the sum of two 
parts Ec — 6flei + SQspin, where 
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Here G^ujm) and F(uJm) are the real and anomalous 
Greens functions given by Gu;^ (k) — — (ek + i^uj„^ ) /(^k + 



+ 'J'D and - *<i>(c^™)/(e^ + S^^, 

where ek = ^^(fc — kp), and Ec^^ = ujm + 



. The 

conventionally defined pairing gap A^^^^ is the ratio 
of the anomalous vertex and the self-energy: A^^^ = 
^w„Wm/Sij^^. In the second term, x (Qi ^) is the dynam- 
ical spin susceptibility related to the spin polarization 
operator H^^ by x (q, ^m) = x{Q, 0)/(l + (q - Q)^ - 
n^,^) where xiQi 0) is the static staggered susceptibility. 
The factor of 3 in Qspin is due to the spin summation. 

The term Silei — f^gf — f^"; accounts explicitly for the 
appearance of the anomalous pairing vertex <i>„^ , and for 
the feedback changes to the fermionic self-energy. The 
term S^spin = ^Ipin ~ ^spm accounts for changes to the 
spin propagator via the changes to the spin polarization 
operator Iluj„^ ■ This term is almost negligible for phonon 
superconductors [ict| . It is however important for mag- 
netic superconductors where the bosonic mode that me- 
diates the superconductivity is by itself affected by the 
pairing. 

The fermionic self-energy , the pairing vertex ^uj^ , 
and the spin polarization operator are the solutions 
of three coupled Eliashberg equations |ll| 
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Eqs. (|^) contain only two inputs: the overall energy 
scale u) that is set by the spin-fermion interaction, and 
the dimensionless spin-fermion coupling A oc ^ that di- 
verges as the system approaches the antiferromagnetic 
instability. Physically, the energy scale uj is the ultimate 
upper cutoff for the strong coupling behavior (S„^ < 
for LOm > Wm), while dimensionless A can be represented 
as the ratio (2A)^ = uj/uJsf of u) and another typical scale 
LUsf that sets the upper boundary of the Fermi-liquid be- 
havior in the normal state. 

The evolution of the condensation energy with A can be 
understood as follows. At small A, i.e., at weak coupling, 
w < iOsf , and the system behaves as a conventional Fermi 
liquid. In this limit, the pairing potential is static, i.e., 
the condensation energy is entirely electronic and BCS 
like. At strong couplings, i.e. A > 1, a) > Ws/, there is 



a frequency range between w^/ and Co within which the 
system behaves as a non Fermi liquid. Previous studies 
demonstrated that in this regime both the onset of 
the pairing instability and the magnitude of A„ at T = 
are determined by u), and hence the condensation energy 
results predominantly from fermions located in the non- 
Fermi liquid frequency range. For these fermions, retar- 
dation effects not included in the BCS theory become 
dominant, hence one expects strong deviations from the 
BCS form of Ec- Comparison to experiment jj, |ll| yields 
A ~ 1.5 — 2 at optimal doping. This doping already cor- 
responds to a strong coupling (iD ~ (10 — 16)ujsf). 

As S and $ depend only on w, the momentum integra- 
tion in Eqs. (pi) can be performed explicitly and yields 
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where 11,, 



In the normal state, li^ 



-\ujm\- We introduce TV/ via / (fk/{'i-K'^) = Nf J dck and 
an analogous spin "density of states" Ng = {8/3n'^) Oj/vp, 
where vp is the Fermi velocity. One can make sure 
that for A > 1, the frequency sums in both terms yield 
0(((D)2), hence Silgi ~ ^f^"^, while SQsptn ~ NsU)^. 

In Fig.|l| we present the results for the electronic and 
spin contributions to the condensation energy for differ- 
ent values of the coupling A. We used previously obtained 
solutions of the Eliashberg equations at T « in real fre- 
quencies, and related 6fl with Ei^, and along the 
real frequency axis using the spectral representation. 

There are three striking features of Figj^. First, the 
electronic contribution to the condensation energy is pos- 
itive. As in the BCS limit the electronic condensation 
energy is negative and equal to —NfA'^/2, this implies 
that the electronic condensation energy changes sign at 
a rather small A ^ 0.4, and is positive for all A > 1/2 
presented in the figure. Second, for all A shown, the spin 
part Sflspin is negative. Third, at large A, the two parts 
of the condensation energy nearly saturate. 

We next have to find a relation between Ng and Nf . 
We use a trick employed by Bardeen and Stephen Q 
and introduce a "mixed" self-energy S„s that has the 
normal state form, but with the spin propagator from the 
superconducting state. We then introduce the functional 
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FIG. 1: Electronic (Sfl^i), spin {5Qspi„), and total (Ec) con- 
densation energy per unit cell at T = for various couplings 
A. The lines are a guide for the eye. The sum of SQ^i and 
SQapin produces a total condensation energy which is negative. 
We used Nf = Ist/eV, and Ns ~ 0.17st/eV as explained in 
the text. 

the fermionic self-energy via the spin polarization oper- 
ator and then integrating over the bosonic momentum. 
The exphcit integration yields 

y = 7rTiV/^(|S„,^„|-|S„,,^„J) (5) 

ni 

Expressing E in terms of 11 , and using the fact that 
the polarization operator made out of Gns is precisely 
n„(w) as the latter depends on the functional form of 
the Green's function but not on the magnitude of the 
self-energy , we find after simple algebra 

Y = -nN, J2 ^n..^ log ""'^'^f--- (6) 

Comparing (||) and we see that Nf and Ng are indeed 
related. We evaluated the ratio for various A and found 
that with small variations, Nf/Ng w 5.9 for all A > 0.5. 

In Fig. |l| we present the result for the full conden- 
sation energy Ec using this ratio of Nf/Ng. To set 
the overall scale, we adopt a commonly used estimate 
Nf = 1st/ eV We emphasize that changing Nf will 
only change the overall scale and not the functional form 
of -Be (A). We see from the figure that the total con- 
densation energy is negative, as it indeed should be in 
a superconductor, but its value is rather small due to a 
substantial cancellation between the spin and electronic 
contributions to Ec- 

We also see that the condensation energy flattens at 
A ~ 2, and decreases at large couplings despite the fact 
that the pairing gap increases monotonically with A [prf . 
The decrease in Ec can be understood as a reflection of 
the fact that as A increases the pairing process more and 
more involves the exchange of classical, on-shell bosons. 
Due to energy conservation, such bosons can not lead 
to a gain in Ec- p3| . This behavior is very counterin- 
tuitive from a BCS perspective, where the condensation 
energy scales with A^. To further emphasize this point. 
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FIG. 2: Total condensation energy Ec compared with the 
BCS result Sf^^ at T = for various couplings A. We 
used Nf = Ist/eV and Nf/Ns ~ 5.9 as explained in the 
text. Observe that the BCS condensation energy increases 
monotonically as the coupling gets larger, while the actual 
condensation energy flattens at A ~ 2 and decreases slightly 
at large couplings. 

we plot in Fig.||our Ec and the BCS condensation energy 
— with the same A and Nf. We clearly see that 
for A > 1, corresponding to optimally doped and under- 
doped cuprates, BCS theory yields qualitatively incorrect 
results for Ec- 

We next present the results for the change in the ki- 
netic energy when the system enters the superconducting 
state. The conventionally defined kinetic- energy for an 
interacting fermionic system is 

E,,^ = 2TJ2 J^e,G^^{k) (7) 

where G^j^ (k) is the full fermionic Green's function that 
contains the self-energy. Integrating over momentum and 
subtracting the normal state result from Ekm in a super- 
conductor we obtain 

6Eun = 2Nf7TTY,\/WZ+^-\^n,.J («) 

m 

In the BCS hmit, = A, S = 0, = w„ and 

SEkin = 2Ar/7rTX;,„ + - \^™ \ which is obvi- 

ously positive. The results at finite A are presented in 
Fig ||. We used the same computational procedure as 
before. At low couplings the kinetic energy is positive, 
as one naively expects. At larger A, however, the kinetic 
energy passes through a maximum at A 2 and then 
becomes negative at large A. This implies, as we an- 
ticipated, that at strong coupling, the lowering of i?fei„ 
via the change in the self-energy due to the "undressing" 
of fermions overcomes the effect of particle-hole mixing 
tliat increases Ekm- This behaviour is very similar to 
that obtained by Norman et. al. 

As the first term in Silei is equal to —SEkin/2 (eqs. 
^ and H), one can indeed argue that the condensation 
energy at large couplings is at least partly driven by the 
lowering of the kinetic energy. However, a comparison 
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FIG. 3: Kinetic energy SE^in compared with total condensa- 
tion energy Ec a,t T — for various couplings A. The param- 
eters are the same as in Fig. 2. The kinetic energy change is 
positive at low couplings, but negative at high coupling. 

of Figs. |l] and || shows that this is just another way 
to interpret strong coupling effects that affect both the 
fermionic and bosonic propagators via mutual feedback. 

We emphasize that the calculated value of Ec is small, 
only ~ 15K around optimal doping. This is rather re- 
markable as all typical energies in the problem are much 
higher, i.e oj - (2.5 - 3) lO^K. This small value of Ec 
is partly due to small prefactors, but is also the result 
of substantial cancellation between the spin and elec- 
tronic contributions to Ec. Loram et al extracted 
Ec ~ 0. 12fcB Tc ^ 10 K from the jump of the specific heat 
at Tc- Our result is is very close to this value. 

Loram et al found that the condensation energy drops 
as one moves into the undcrdoped regime We also 
found a decrease in Ec at large couplings, although not as 
spectacular as in Loram's experiments. This is due to the 
strong dispersion of the bosonic propagator assumed in 
the spin-fermion model. A weaker spin excitation disper- 
sion, as suggested by neutron experiments fl^ , should 
result in a stronger reduction of Ec- 

Finally, we discuss the accuracy of the Eliashberg ap- 
proach and also clarify what we mean by strong coupling. 
As discussed previously the validity of Eliashberg 



theory is related to the development of an "effective" 
Migdal theorem at strong coupling allowing vertex cor- 
rections to be neglected. As spin fluctuations become 
overdamped due to spin-fermion exchange they become 
slow compared to fermions (just as phonons are slow 
compared to electrons). This effective Migdal theorem 
emerges if the momentum integration is confined to the 
vicinity of the hot spots which in turn requires Co to be 
smaller than the fermionic bandwidth W. This defines 
a "universal" strong coupling limit: X > l,uj < W. The 
neglect of the A:— dependence of the d— wave gap and of 
the self-energy along the Fermi surface (i.e., the approx- 
imation of Nf by a constant) is an 0(1) approximation 
even when ui ^W, but these two /c— dependences are to- 
tally irrelevant from the physics perspective and can be 
safely omitted |^]. Furthermore, even when u) > 
the pairing problem does not change qualitatively ex- 
cept that typical momenta measured from a hot spot 
are now 0{pf), and and Nf are no longer universally 
related. In this regime, the spin part of the condensation 
energy prevails, and scales as W'^/ui ~ J, as predicted by 
Scalapino and White |3|. 

To conclude, we have shown that at strong spin- 
fermion coupling, the behavior of the condensation en- 
ergy differs strongly from the BCS prediction. Instead of 
getting larger as A increases, as predicted by BCS, the 
total condensation energy flattens off and then decreases. 
This decrease appears to be correlated to a change in 
sign of the kinetic energy, another non-BCS type behav- 
ior. We also found that the spin and charge contributions 
to Ec are of comparable strength and opposite sign, and 
negative Ec results from a delicate balance between the 
two. This behavior has no analog for phonon supercon- 
ductors, where the feedback on bosons is a minor effect. 
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